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I.  Program  goal 


The  goals  for  the  project  are  to  solve  matrix  algebraic  problems  at  high  speed  by 
applying  photorefractive  crystals  in  optical  processing  systems,  and  to  study  the 
photorefractives  for  such  application.  Although  our  interest  was  primarily  on  solving 
matrix  algebraic  problems,  in  order  to  improve  system  performance  we  found  it  important 
to  optimize  crystal  operation  and  to  improve  the  photorefractive  memory  performance. 
For  example,  a  fast  response  and  large  coupling  gain  bandwidth  photorefractive  amplifier 
is  required  in  the  matrix  inversion  system,  and  high  photorefractive  efficiency  and  large 
storage  capacity  are  required  for  the  correlation  matrix-tensor  multiplier  system. 

EL  Program  Summary 

The  report  will  be  divided  into  four  sections.  Section  2.1  describes  an  all-optical 
implementation  of  an  iterative  algorithm  for  matrix  inversion.  Section  2.2  reports  the 
implementation  of  a  correlation  matrix  tensor  multiplier  algorithm  in  photorefractive 
crystals.  Section  2.3  summaries  the  results  that  have  been  obtained  in  improving 
photorefractive  performance,  which  includes  45°-cut  BaTiC>3,  applying  electrical  field 
and  moving  grating  on  SBN  and  GaP.  Section  2.4  reviews  our  work  on  photorefractive 
memory,  which  includes  circulating  memory,  incremental  recording,  orthogonal  phase 
encoding  and  selective  erasure.  A  list  of  publications  resulted  from  the  research  is 
included  at  the  end  of  this  report. 


2.1  Optical  implementation  of  an  iterative  algorithm  for  matrix  inversion 

Many  areas  of  engineering  require  real-time  matrix  inversion.  It  is  often  found  that 
matrix  inversion  are  time-consuming  for  digital  electronic  processors.  Therefore,  parallel 
optical  processing  approaches  to  the  solution  are  of  great  interest  Because  the  solution 
require  iterative  algorithms  which  is  difficult  to  achieve  by  electronic  processors  for  large 
size  matrices,  the  use  of  feedback  techniques  in  optical  information  processing  is 
extremely  important  We  have  carried  out  an  all-optical  implementations  of  an  iterative 
algorithm  for  matrix  inversion  (Ref.l)  using  a  confocal  Fabry-Perot  system  with  optical 
feedback  and  gain  provided  by  a  photorefractive  BaTiC>3  crystal  in  the  feedback  path. 

Fig.l  shows  the  schematic  diagram  of  the  matrix  inversion  processor.  Y  and  X 
represent  input  and  output  vectors  of  the  processor,  respectively.  In  the  feedback  loop,  a 
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matrix- vector  multiplier  forms  the  product  AX.  The  input  and  output  of  the  processor  can 
be  related  according  to  the  equation 
Y±AX  =  X. 

The  output  vector  of  the  processor  is  therefore,  given  by 
x-  (i±A)-iy  =b -*y. 

To  find  the  inverse  of  the  matrix  B,  we  incorporate  the  corresponding  matrix  A  =  I  ;  B  in 
the  optical  processor.  When  the  input  vector  Y  is  set  to  zero  everywhere  except  for  one 
element,  whose  value  is  normalized  to  unity,  and  the  position  of  this  unity  element  in  Y  is 
sequentially  shifted  from  the  uppermost  to  the  lowermost  position,  the  matrix  B*1  =  (I  - 
A)-1  will  be  generated  in  real-time  at  the  output,  column  by  column. 


Fig.l  Schematic  diagram  of  the  matrix  inversion  processor  showing 
the  matrix-vector  multiplier  in  the  feedback  loop.  Y  and  X  represent 
(Nxl)  input  and  output  vectors  respectively. 
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photorefractive  image  amplifier  into  the  feedback  loop  of  the  confocal 
Fabry-perot  processor  made  of  mirrors  Mi  and  M2;  L  and  C  are  spherical 
and  cylindrical  lenses  respectively,  PR  -  BaTi03  Photorefractive  crystal 
D-  Dove  prism,  S  -  slit 


Fig.2  shows  the  optical  system  corresponding  to  the  diagram  of  Fig.l  for 
implementing  the  matrix  inversion  algorithm,  with  the  matrix-vector  multiplier  and  the 
photorefractive  amplifier  incorporated  into  a  confocal  Fabry-Perot  processor.  The  input 
slit  S0  which  represents  the  nonzero  element  of  the  input  vector  Y  was  rotated  by  the 
Dove  prism  D  by  90°  and  was  then  telecentrically  imaged  to  the  matrix  mask  A  by  lenses 
Lq-Lq’  and  spherical  mirror  Mi  and  M2.  The  slit  S  selects  the  vertical  dc  component  of 
the  Fourier  Transform  of  AX..  A  photorefractive  crystal  of  regular  cut  BaTiC>3  (crystal 
faces  perpendicular  to  crystal  c-axis)  was  introduced  in  the  feedback  loop  as  a  gain 
media.  The  output  B*1  was  observed  through  output  lenses  L4C4  column  by  column  as 
the  input  slit  SQ  is  spatially  shifted  sequentially. 


Fig.3.  Photograph  of  the  experimental  setup.  The  spatial  filters 
(matrix  A  and  slit  S)  are  inserted  into  the  liquid  gate.  The  minor 
on  the  right  is  mounted  on  a  piezoelectric  translator. 

Fig.3  shows  the  photograph  of  the  experimental  setup.  One  of  the  two  spherical 
minors  is  mounted  on  a  piezoelectric  translator  which  permits  fine  adjustment  of  the 
mirror  spacing  allowing  the  selection  of  either  a  positive  or  negative  feedback 
configuration.  A  2-D  gray  level  mask  representing  matrix  A  of  size  3x3-mm  is  generated 
on  spectroscopic  649F  Kodak  film  via  a  computer-controlled  scanner.  The  slit  S  of  width 
10  pm  is  used  to  ensure  the  maximum  relative  error  for  the  matrix- vector  multiplier  to  be 
<  5%.  A  single  crystal  of  BaTi03  was  used  for  coherent  image  amplification.  This 
amplifier  was  optimized  with  respect  to  such  operational  characteristics  as  gain  vs  space- 


bandwidth  product  and  dynamic  range.  The  gain  of  the  photorefractive  amplifier  is  also 
required  to  be  constant  over  the  entire  dynamic  range  of  the  matrix  A.  This  condition 
was  fulfilled  by  setting  the  pump  beam  intensity  much  larger  than  the  signal  beam 
intensity.  A  gray  level  dynamic  range  of  the  order  of  100-500  was  demonstrated  for 
BaTiC>3  crystals.  A  large  space-bandwidth  product  of  106  was  also  demonstrated; 
therefore  the  system  has  the  potential  to  invert  very  large  size  matrices  (103  x  103). 
Finally,  the  processing  time  is  independent  of  the  size  of  the  matrix  and  is  only 
determined  by  the  time  response  of  the  crystal,  which  is  inversely  proportional  to  the  total 
incident  intensity  Iq,  typically  0.1-second  to  1-second  for  Io  =  100  to  10  mW-cnr2. 

Experimental  verification  was  performed  using  an  (8  x  8)  matrix.  Fig.4  shows  the 
experimental  result  provided  by  a  3-D  plot  of  the  intensity  of  the  output  matrix  B'1. 
Comparing  with  the  digital  computation  result  of 
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the  experimental  accuracy  is  5%. 


•**-»  s: 


Hg.4.  Experimental  result  of  the  inverted  matrix  B'1  represented 
in  three-dimension. 
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2.2  Optical  Implementation  of  Correlation  Matrix-Tensor  Multiplier  Algorithm 

Matrix-tensor  multiplication  has  many  important  applications,  including 
interconnection  between  2-D  data  arrays  for  digital  and  neural  computing.  For  NxN  input 
and  output  arrays,  the  interconnection  tensor  is  4th  rank  containing  the  N4  independent 
values . 

We  have  built  a  correlation  matrix-tensor  multiplier  (CMTM)  that  can  implement  an 
arbitrary  complex-weighted  2-D  array  interconnection.  In  the  CMTM  system,  the  2-D 
input  was  transmitted  through  a  phase-code  mask  and  correlated  with  a  control  image 
constructed  from  the  interconnection  weights.  The  optical  correlation  was  performed  by 
either  photorefractive  real-time  four-wave  mixing  or  by  pre-storing  the  interconnection 
weight  in  the  crystal.  Because  the  real-time  four-wave  mixing  is  limited  by  slow 
photorefractive  response  (Ref.3),  we  chose  to  use  the  late  version  for  implementation. 
Multiple  interconnection  patterns  were  prestored  as  color-multiplexed  volume  reflection 
holograms  in  z-cut  LiNbC>3.  Fast  reconfiguration  between  interconnection  patterns  is 
possible  using  a  wavelength  tunable  source,  decoupling  both  the  modulation  and 
switching  speeds  from  slow  photorefractive  response.  The  CMTM  algorithm  allows  for 
graceful  accommodation  of  limited  spatial  light  modulator  size  needed  to  input  the 
control  tensor  by  trading-off  output  signal-to-noise  ratio  for  reduced  control  tensor 
bandwidth  requirement. (Ref . 2)  Experimental  results  confirmed  theoretical  predictions 
that  the  CMTM  system  works  best  for  densely  connected  networks,  with  a  large  fan-in  to 
each  output  Interconnection  of  up  to  4096  inputs  and  outputs  was  demonstrated  using 
such  dense  interconnection  patterns.  An  aggregate  average  SNR  of  over  200  was 
obtained  for  1024  inputs  and  outputs.  Fig.5  shows  the  output  for  different  input  array 
size  (N2)  and  phase-code  density  k.  The  connection  tensor  size  was  1000x1000  pixels 
(25-mm  square),  approaching  the  realistic  maximum  for  an  optical  system  to  resolve  with 
good  accuracy.  The  overall  results  match  the  theoretical  prediction. 
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Fig.5.  Effect  of  array  size  on  output  SNR.  The  control  image  SBP 
was  held  at  about  1000x1000  by  decreasing  k  linearly  with  N.  At  left, 
all  the  inputs  are  on.  At  right,  half  the  inputs  are  masked.  Output  SNR 
decreased  as  N2  grows  from  256  to  1024  to  4096. 
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2.3  Techniques  developed  at  UCSD  for  improving  photorefractive  performance 
The  system  performance  of  both  matrix  inversion  and  CMTM  strongly  relay  upon  the 
photorefractive  operation.  For  example,  a  fast  response  and  large  coupling  gain 
bandwidth  photorefractive  amplifier  is  required  in  the  matrix  inversion  system,  and  high 
photorefractive  diffraction  efficiency  and  large  storage  capacity  are  required  in  the 
CMTM  system.  In  the  following  we  summarize  some  of  the  techniques  of 
photorefractive  performance  improvement  developed  at  UCSD  under  the  support  of  this 
contract 

2.3.1  Enhanced  photorefractive  performance  from  45°-cut  BaTiC>3 

BaTiC>3  photorefractive  amplifier  is  one  of  the  key  elements  in  the  matrix  inversion 
implementation  reported  in  section  2.1.  In  this  implementation,  the  BaTi03  is  a 
conventional  cut  with  its  faces  polished  perpendicular  to  the  crystal  optic  axis  c .  In  order 
to  improve  the  system  performance  such  as  gain,  response  time  and  uniformity,  we 
proposed  and  have  demonstrated  that  changing  the  crystal  cut  of  BaTiC>3  to  45°  to  c  can 
provide  superior  photorefractive  performance  to  our  system  applications.(Ref.3)  It 
increases  the  achievable  gain  coefficient,  reduces  the  response  time,  and  changes  the  gain 
and  index  modulation  response  as  a  function  of  external  input  angles. 

We  characterized  the  photorefractive  performance  of  a  45°-cut  BaTiC>3  and  compared 
it  with  the  conventional  0°-cut  material.  The  achievable  peak  gain  coefficient  T  is 
increased  to  a  value  of  26  cm-1,  while  T  for  a  conventional  BaTiC>3  is  only  12  cm*1.  The 
strong  fanning  observed  in  thicker  45°  crystal  reduces  the  effective  T  for  two-wave 
mixing  but  can  be  applied  to  enhance  the  performance  of  the  crystal  in  certain 
photorefractive  applications  which  use  the  fanning  effect  to  their  advantage.  Two 
examples  of  these  are  self-pumped  phase  conjugation  and  fanning-based  novelty  filtering 
(Ref.4).  The  angular  dependence  of  T  on  the  input  beam  angles  is  flatter  in  the  45°-cut 
crystal,  resulting  in  more  uniform  amplification  over  a  wide  range  of  input  beam  angles. 
The  response  times  measured  in  45°-cut  crystal  were  1  order  of  magnitude  faster  than 
similar  0°-cut  crystals.  Application  of  an  external  electric  field  can  slightly  increase  I\ 
but  does  not  reduce  the  response  time  and  can  damage  the  crystal. 

2.3 X  Increased  multiplex  hologram  capacity  of  SBN  with  applied  field 

Increasing  storage  capacity  of  photorefractive  materials  is  very  important  for 
improving  the  CMTM  performance  because  more  interconnection  patterns  can  be 
prestored  and  accessed  at  very  high  speed.  We  have  characterized  the  performance  of 
photorefractive  Sro.6Bao.4Nb206  under  an  externally  applied  electric  field.  We 
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demonstrated  that  the  volume  multiplexing  capacity  increases  substantially  by  using  an 
applied  dc  field.(Ref.5)  In  a  1-mm-thick  crystal  the  number  of  2%  diffraction  efficiency 
holograms  increased  from  4  to  32  at  an  external  field  of  14  kV/cm.  The  minimum 
diffraction  efficiency  for  30  holograms  increased  from  2.4%  to  5.4%  at  12.6kV/cm  field. 
Accurate  prediction  of  photorefractive  performance  under  an  applied  field  can  be  made 
with  material  parameters  measured  at  zero  applied  field.  The  field  dependence  of  the 
gain  coefficient,  response  time,  sensitivity,  and  diffraction  efficiency  were  obtained  and 
fit  very  well  with  the  theoretical  prediction  calculated  using  Kukhtarev’s  photorefractive 
material  equation.  Optimum  crystal  operation  conditions  of  applied  field  and  grating 
period  were  obtained. 

2.3  J.  Moving  grating  for  enhanced  holographic  recording  in  SBN 

Energy  coupling  during  holographic  recording  in  thick  SBN  (e.g.  5-mm)  tends  to 
reduce  the  beam  ratio  of  1  at  the  input  to  10" 7  at  the  exit  face.  This  greatly  limits  the 
usable  crystal  thickness  reducing  the  diffraction  efficiency  and  storage  capacity  which 
are  required  by  many  applications,  including  CMTM.  We  have  demonstrated  that 
holographic  recording  can  be  improved  by  using  a  moving  grating  to  reduce  energy 
coupling.(Ref.6)  In  Ce*doped  Sro.eBao.4Nb206  under  an  applied  dc  electric  field, 
experimental  and  theoretical  results  show  that  the  proper  grating  velocity  can  maximize 
the  modulus  and  the  real  part  of  the  space-charge  field  while  reducing  the  imaginary 
component  of  the  space-charge  field  (and  energy  coupling)  to  zero.  Avoiding  energy 
transfer  between  the  recording  beams  allows  us  to  maintain  maximum  contrast 
throughout  the  crystal,  producing  a  uniform,  high-index-modulation  grating  with 
enhanced  diffraction  efficiency  and  superposition  properties.  We  demonstrated  that  with 
an  applied  field  of  5.3kV/cm,  the  peak  diffraction  efficiency  was  increased  from  58%  to 
93%  by  choosing  an  optimum  moving  grating  velocity  which  minimized  coupling. 

2.3.4.  Moving  grating  for  enhanced  beam-coupling  and  phase-conjugation  in  GaP 
The  ferroelectric  photorefractive  materials,  such  as  BaTi03,  SBN  and  LiNbC>3 
possess  slow  time  response,  affecting  the  speed  of  both  matrix  inversion  system  and  real¬ 
time  CMTM  system.  In  order  to  improve  the  system  operating  speed,  we  explored  new 
photorefractive  materials  with  fast  time  response.  It  has  been  shown  that  semiconductor 
photorefractive  materials  provide  response  time  much  fast  than  that  of  ferroelectric 
materials,  among  which  GaP  is  particularly  attractive  because  it  operates  in  the  spectral 
range  of  0.6  to  0.9  pm.  However,  weak  photorefractive  effect  in  GaP  has  been  limiting 
its  applications.  We  have  demonstrated  a  technique  to  enhanced  photorefractive  effect  in 
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GaP  by  using  a  dc  field  and  a  moving  grating  at  laser  wavelength  of  633nm.  Two-beam 
coupling  gain  coefficient  T  was  increased  from  0.3  cm*1  to  2.5  cm*1  and  the  transient 
phase-conjugate  reflectivity  was  increased  from  0.3%  to  4.5%  by  applying  a  dc  field  of 
16.7  kV/cm  and  choosing  the  optimum  fringe  velocity.  The  experimental  results  are 
found  in  good  agreement  with  the  theoretical  predictions.(Ref.7) 

2.4  Photorefractive  Memory 

3-D  photorefractive  memories  offer  the  advantages  of  high  storage  density, 
distributed  storage,  and  fast  parallel  access.  These  characteristics  have  been  applied  to 
our  CMTM  system  (Section  2.2)  and  many  other  parallel  opto-clectronic  processing 
applications.  The  high  performance  of  3-D  photorefractive  memories  rely  on  various 
basic  requirements,  such  as  fast  random  access  time,  read/write/erase  capabilities,  high 
angular  and  wavelength  selectivity,  high  resolution,  large  and  uniform  diffraction 
efficiency,  long  storage  time  and  nondestructive  readout,  which  are  still  the  research 
topics  and  are  extremely  important  to  our  CMTM  system  performance.  In  this  section, 
we  will  review  the  related  research  projects  which  have  been  accomplished  at  UCSD 
during  the  last  few  years. 

2.4.1.  Dynamic  circulating  photorefractive  memory 

The  objective  of  this  project  is  to  construct  a  dynamic  photorefractive  memory  that 
provides  (1)  nondestructive  readout,  i.e.,  after  many  readings,  the  memory  can  be 
refreshed  by  rewriting  the  still-useful  information  from  one  crystal  into  the  other,  and  (2) 
write-erase  capability,  i.e.,  while  rewriting  the  useful  part  of  the  information  from  one 
crystal  to  the  other,  we  also  can  replace  obsolete  information  that  must  be  erased  by  new 
information  introduced  from  the  memory  input  For  accomplishing  this  task,  we  proposed 
and  have  implemented  an  architectural  approach  using  dynamic  3-D  volume  storage  that 
circulates  holograms  between  two  photorefractive  crystals.(Ref.8) 

In  our  architecture  the  input  images  are  stored  sequentially  in  one  of  the  two  crystals 
(i.e.,  PRC1)  until  the  total  memory  capacity  of  that  crystal  is  reached.  After  many 
readouts  the  stored  information  in  PRC1  will  be  degraded.  We  then  transfer  the  useful 
part  of  the  information  from  PRC1  to  the  second  crystal  PRC2.  If  the  storage  capacity  of 
PRC2  is  not  reached  after  transfer,  additional  input  images  can  be  recorded  in  PRC2.  The 
first  crystal  is  then  optically  erased  to  prepare  it  for  the  next  transfer  cycle.  The  two 
crystals  are  functionally  identical,  exchanging  tasks  periodically.  In  order  to  increase  the 
effective  write-erase  asymmetry  of  the  crystal  and  amplify  the  recalled  images,  an  optical 
amplifier  (LCLV)  was  introduced  into  the  system.  We  demonstrated  transfer  of  a  single 
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image  (piwel  size  400  jxm)  between  two  crystals.  No  noticeable  degradation  in  the  pixel 
qual:v  was  observed  after  more  than  fifteen  transfers  with  the  average  SNR  better  than 
86.  This  experiment  shown  that  our  circulating  memory  architecture  provides  write-erase 
storage  capacity  with  robustness  to  multiple  optical  readouts. 

2.4.2.  Incremental  recording  and  orthogonal  phase  codes  for  photorefractive 
hologram  multiplexing 

The  commonly  used  schedule  recording  technique  for  photorefractive  hologram 
superposition  involves  a  complicated  recording  procedure  with  a  long  first  exposure 
followed  by  shorter  and  shorter  exposures.  The  schedule  is  calculated  from  the  material’s 
response  times  and  the  maximum  attainable  index  modulation.  A  small  error  in  material 
characterization  (also  nonuniformities  or  changes  in  material  characteristics)  can  result  in 
highly  nonuniform  diffraction  efficiencies.  Also  problems  of  fanning  and  beam  coupling 
will  be  dominant  in  the  long  initial  exposures,  which  tend  to  restrict  the  maximum 
attainable  index  modulation  to  well  below  the  theoretical  value.  We  have  demonstrated  a 
new  recording  schedule  called  incremental  recording  by  which  the  problems  associated 
with  the  schedule  recording  can  be  avoided.  We  have  studied  the  effecis  of  the 
incremental  recording  approach  on  holographic  storage,  and  the  criteria  for  choosing 
appropriate  recording  parameters.  In  this  approach,  the  diffraction  efficiency  uniformity 
dependents  only  on  At  but  not  on  material  parameters.  Therefore,  by  using  incremental 
recording  it  is  much  easier  to  reach  uniform  diffraction  efficiency  than  by  using  schedule 
recording  (Ref.9). 

In  incremental  recording,  each  of  N  holograms  is  recorded  with  a  series  of  exposures 
At,  each  extremely  short  compared  with  the  material’s  response  time.  During  recording, 
each  image  and  reference  pair  is  sequentially  displayed,  repetitively  cycling  through  all  N 
images.  The  first  recorded  hologram  will  remain  after  incrementing  the  other  (N-l) 
holograms  because  small  index  modulation  possess  high  writing  and  low  erasing 
sensitivity.  As  this  process  is  repeated,  all  the  holograms’  diffraction  efficiencies 
gradually  increase  until  the  recording  process  reaches  saturation,  i.e.,  when  index  gained 
during  recording  increment  equals  to  that  erased  during  (N-l)  erasures  when  the 
holograms  are  being  recorded. 

In  incremental  recording  approach,  the  recording  system  must  repeatedly  display  the 
input  and  reference  image  pairs,  reproducing  the  interference  pattern  to  within  much  less 
than  one  fringe  width.  To  achieve  this  we  have  developed  a  method  of  addressing  called 
orthogonal  phase  encoding  by  which  different  holograms  can  be  recorded  and  addressed 
without  mechanical  movement  In  orthogonal  phase  encoding  approach  the  reference 
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wave  fronts  are  generated  by  displaying  a  set  of  orthogonal  phase  codes  on  a  phase- 
modulator  array.  Each  hologram  is  recorded  (addressed)  by  one  of  the  set  of  phase  coded 
plane  waves.  The  phase  codes  are  chosen  such  that  the  inner  product  of  one  code  with 
any  other  code  from  the  set  is  zero,  that  is  ,  they  are  orthogonal  to  each  other.  When  the 
crystal  is  illuminated  with  one  phase  code,  the  unwanted  cross  terms  from  the  other 
holograms  interfere  destructively  in  the  volume  of  the  crystal,  while  the  output  recorded 
with  that  phase  code  interference  constructively  to  produce  a  full-intensity  reconstruction. 
In  general,  a  phase-modulator  with  M  pixels  provides  M  orthogonal  phase  codes.  In  the 
particular  case  of  M  =  2k,  the  solution  is  a  set  of  binary  phase  codes  equivalent  to  the 
bipolar  Walsh  function.  These  codes  can  be  generated  from  a  k  =  log2M  bit  address, 
therefore  this  phase-encoding  technique  is  useful  for  fast  parallel  access,  especially  when 
a  smart  spatial  light  phase  modulator  with  local  electronic  logic  is  used. 

The  combination  of  incremental  recording  and  binary  orthogonal  phase  codes  was 
demonstrated  experimentally.  We  tested  the  orthogonality  of  codes  by  incrementally 
recording  16  images  in  a  1-mm-thick  Fe-doped  LiNb03  crystal.  The  cross  talk  between 
the  images  was  low,  providing  an  average  SNR  of  100:1.  We  also  investigated  the  signal 
uniformity  for  the  incremental  recording  by  varying  the  incremental  exposure  time.  As 
expected  by  our  theoretical  prediction,  the  signal  uniformity  is  significantly  improved 
with  a  reduced  recording  increment  The  storage  capacity  of  the  incremental  recording 
procedure  is  identical  to  that  of  the  schedule  recording. 

2.43.  Selective  erasure  for  dynamic  photorefractive  memory 

Applications  of  photorefractive  memory  require  not  only  robust  readout  capability 
and  uniform  diffraction  efficiency  with  high  signal  Fidelity  which  are  the  subject  of 
Section  2.4.1  and  Section  2.4.2  respectively,  but  also  fast  update  capability  for  such 
applications  as  neural  networks.  The  simple  way  of  incoherent  erasure  which  has  been 
commonly  used  is  too  slow,  while  selective  erasure  based  on  superimposed  recording  of  a 
n-phase-shifted  image  onto  the  old  image  has  been  reported  to  be  much  faster  than  the 

incoherent  erasure.  The  objective  of  this  project  is  to  analyze  the  dynamics  of  the 
selective  erasure  process  and  to  apply  it  to  a  dynamic  photorefractive  memory. 

We  derived  a  theoretical  equation  of  the  dynamics  of  the  index  distribution  for 
selective  erasure  process  by  using  photorefractive  material  equations.  This  equation  can 
be  used  to  predict  the  diffraction  efficiency  versus  erasure  time  and  selective  erasure  time 
(the  time  for  completely  erasing  the  hologram)  versus  initial  refractive  index  modulation. 
It  was  shown  that  selective  erasure  process  is  much  faster  than  incoherent  erasure 
process.  Experimental  results  show  good  agreement  with  theoretical  prediction.  (Ref.  10) 
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We  have  applied  the  selective  erasure  process  together  with  the  incremental  recording 
schedule  (Section  2.4.2)  to  the  dynamic  photorefractive  optical  memory.  Four  pages  of 
information  were  stored  in  a  1-mm-thick  Ce-doped  Sro.6Bao.4Nb206  by  incremental 
recording  to  reach  saturation.  Fig.6  shows  the  experimental  demonstration  of  selective 
erasure  process.  One  of  the  four  recorded  binary  images  of  the  letter  ( U,  C,  S,  D  )  [see 
Fig.6(a)]  D  was  selectively  erased  [see  Fig.6(b)]  and  was  replaced  by  a  C  [see  Fig.6(c)]. 
It  took  only  a  few  seconds  to  completely  selective  erase  one  page  and  another  20-30 
seconds  to  recording  a  new  page  and  circulating  all  the  pages  until  reach  saturated 
diffraction  efficiency.  The  experiments  iemonstrated  the  robust  readout  and  fast 
read/write/erase  capability  (Ref.  11, 12). 


Fig.6.  Photographs  of  stored  binary  images  (a)  before  selective  erasure, 
(b)  after  selective  erasure,  (c)  after  image  update. 


m.  Conclusion  and  Future  Works 

We  have  demonstrated  an  all-optical  implementation  of  an  iterative  algorithm  for 
matrix  inversion  and  32x32  interconnection  by  using  correlation  matrix-tensor  multiplier 
algorithm.  We  made  several  contributions  to  photorefractive  performance  improvement, 
which  includes  enhanced  gain  and  response  time  using  45°-cut  BaTiC>3,  increased 
holographic  capacity  by  applying  dc  field  on  SBN,  increased  diffraction  efficiency  in 
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SBN  by  applying  moving  grating,  and  enhanced  gain  and  phase-conjugation  in  GaP  by 
applying  dc  field  and  moving  gratirg.  For  photorefractive  memory,  we  have  introduced 
solutions  to  some  of  the  most  serious  problems  such  as  nondestructive  readout, 
write/erase/rewrite  capability,  uniformity  of  diffraction  efficiency  and  fast  random  access 
by  demonstrating  dynamic  circulating  memory,  incremental  recording,  orthogonal  phase 
encoding  address  and  selective  erasure. 

Our  future  work  will  be  more  concentrated  on  applying  photorefractive  materials  for 
packaging  of  opto-electronic  system  such  as  neural  networks  and  optical  interconnection 
systems  in  order  to  fully  use  the  in-situ  recording  nature  of  the  photoreffactives  for  self¬ 
alignment  and  adaptability. 
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